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Применение нечетких чисел для обоснования 
кластеров в методах нечеткой кластеризации” 


В статье рассматривается метод определения оптимального числа кластеров в нечетком с-разбиении, 
основанный на построении интервала значений наиболее возможного числа классов с использованием 
нечетких чисел. Предложена модификация ЕСМ-СУ-алгоритма и приводится результат вычислительного 
эксперимента. 


Введение 


Теория нечетких множеств, предложенная в 1965 году Л.А. Заде [1], нашла 
свое применение практически в любой области научных исследований, в том числе и 
при разработке новых методов автоматической классификации, именуемой также 
кластерным анализом, численной таксономией или распознаванием образов с само- 
обучением [2], в связи с чем Л.А. Заде отмечал, что «глубинная связь между теорией 
нечетких множеств и распознаванием образов основана на том обстоятельстве, что 
большинство реальных классов размыты по своей природе в том смысле, что 
переход от принадлежности к непринадлежности для этих классов скорее постепенен, 
чем скачкообразен» [3]. В работах [4-6] подробно рассматриваются предложенные 
различными исследователями методы нечеткой кластеризации эвристического, 
оптимизационного и иерархического направлений. Оптимизационные методы нечеткого 
подхода к решению задач кластеризации являются наиболее распространенными, 
однако при обращении к указанным кластер-процедурам возникает проблема обоснова- 
ния числа классов, для решения которой традиционно используются различные 
показатели оптимальности числа нечетких кластеров. Вместе с тем, при большом 
количестве предположений о числе классов использование этих показателей также 
затруднительно, и задача состоит в построении множества наиболее возможного 
числа нечетких кластеров. 

Целью предпринятого исследования является решение поставленной задачи, в 
основе которого лежит использование аппарата нечетких чисел. В работе предла- 
гается метод построения допустимого множества значений наиболее возможного числа 
нечетких кластеров в искомом нечетком с-разбиении, основанный на построении 
нечетких величин, исходя из экспертных оценок или результатов разведочного анализа 
данных, а также рассматривается соответствующая модификация предложенного в [7] 
ЕСМ-СУ-алгоритма. Приводятся результаты вычислительного эксперимента и форму- 
лируются предварительные выводы об эффективности предложенного подхода. 


“Исследования проводились при поддержке гранта Президиума Национальной академии наук 
Беларуси в соответствии с Постановлением № 157 Бюро Президиума НАН Беларуси от 11.04.2007. 
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Оптимизационное направление решения нечеткой 
модификации задачи кластеризации 


В нечетких методах автоматической классификации, объединяемых в оптими- 
зационное направление, нечеткая кластеризация понимается как разбиение классифи- 
цируемой совокупности объектов Х = {х.,...,х,„} на семейство его нечетких множеств, так 


что в качестве входного параметра в существующих нечетких методах оптимизационного 
направления автоматической классификации задается число нечетких кластеров с, 
причем при данном подходе под нечетким кластером может пониматься любое нечеткое 
множество, определенное на универсуме. Нечеткие множества ый ‚[=1,....с с соответст- 
вующими функциями принадлежности 4,,..., и, каждого объекта х, , определенные на 
универсуме Х ={х.....х,}, образуют нечеткое с -разбиение, иногда называемое также 
нечетким разбиением в смысле Э.Г. Распини [3], если для каждого объекта х, е Х выпол- 
няется условие у и, =1, и нечеткая модификация задачи автоматической классифи- 
1=1 

кации в экстремальной постановке заключается в нахождении экстремума некоторого 
функционала О(Р) на множестве П всех возможных нечетких с -разбиений Р клас- 
сифицируемого множества объектов Х , что описывается формулой О(Р) > ег. 


Если исходные данные представлены матрицей вида «объект-объект» 
4 „и =[(х,,х,)]=Т...,п,]=1,....п, элементами которой являются попарные коэффи- 
циенты различия между объектами, то для классификации объектов исследуемой 
совокупности используются так называемые реляционные процедуры, основанные 
на минимизации соответствующих функционалов. Примером такого функционала 
может послужить критерий М. Рубенса [8] 
с п Ш 
Оь(Р)=>>. Хирш (к х,}, (1) 
ЕТ 
процедура минимизации которого именуется в литературе ММО?-алгоритмом [8] 
или ЕММ-алгоритмом [6]. В выражении (1) символом 4(х,,х,) обозначается коэф- 


фициент различия между объектами х, их „» Ь]=Ъ....п исследуемой совокупности Х, 


сата(Х) =п,а д, — значение принадлежности 1-го объекта /-му нечеткому кластеру. 
В случае же представления исходных данных в виде матрицы «объект — свойство» 
^^ — [2 ]2=1,....и,Ё =1,....т, задача классификации заключается в минимизации 


критерия, в который введена некоторая метрика, и типичным примером подобных 
функционалов является критерий, предложенный Дж. Данном [9] и позже обобщен- 
ный Дж. Беждеком [10] 


ОРЛ) =ХУик-7 |, (2) 


1=Е = 


процедура минимизации которого широко известна под обозначением ЕСМ-алгоритма. 
Критерий (2), где символом Г’ обозначен прототип /-го нечеткого кластера, а символом 
у — задаваемый исследователем коэффициент нечеткости классификации, такой, что 


1<у<®, послужил основой для целого семейства функционалов и соответствующих 
им нечетких кластер-процедур, которые подробно рассматриваются в [4]. 
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В силу того, что с является параметром любой оптимизационной кластер-про- 
цедуры, одной из главных проблем при использовании оптимизационных методов 
является определение «реального» числа с нечетких кластеров, на которые «рас- 
слаивается» исследуемая совокупность, или, иными словами, проблема обоснования 
числа кластеров, встающая наиболее остро, когда исследователю число классов с 
вообще неизвестно. Для решения этой проблемы были предложены различные пока- 
затели, характеризующие получаемое при использовании того или иного алгоритма 
нечеткое с -разбиение Р”={А',..., А°}.. В частности, для ЕСМ-алгоритма и его модифика- 
ций различными исследователями был введен ряд показателей, наиболее известными 
из которых являются 
— коэффициент разбиения 

1 [Я Ш 5 
7, (Р) => ши, (3) 
П 1 = 
предложенный Дж. Данном [11] и для которого решение задачи определения оптималь- 
ного числа классов в Р” отыскивается в виде тах(7,, (Р)}, с=2,...п-1; 
с 


— энтропия разбиения 


1 Е п 
(р) = ">. ши "п ) (4) 
[Е = 


предложенная Дж. Беждеком, М.П. Уиндхемом и Р. Эрлихом [12], так что оптималь- 
ному числу классов в Р* соответствует ши, (Р)}, с=2,....п-1; 
[6 


— индекс разделимости 


ТУ? 
И ь 
шит 


2 
1 
ж-#| 


(5) 


предложенный Х.Л. Хии Ж. Бени [13], где оптимальное число классов в Р” опреде- 
ляется, исходя из условия тши(И,, (Р)),с=2,....п-1. Для других нечетких кластер- 
[6 


процедур оптимизационного направления также предлагается ряд показателей опти- 
мальности числа классов в нечетком с -разбиении - в частности, в [14] рассматриваются 
показатели оптимальности нечетких кластеров при разбиении исследуемой совокуп- 
ности с помощью ЕММ-алгоритма. Для всех показателей числа классов в нечетком 
с-разбиении решение задачи определения оптимального числа классов в искомом 
нечетком разбиении определяется общим выражением 

ех(У.(Р)), с=2,....п-1, (6) 

с 


где символом Г.(Р) обозначен какой-либо показатель. Таким образом, необходимым 
является проведение серии экспериментов при различных значениях числа классов с, 
для чего оказывается необходимым построение множества С = {с.,...,с”} наиболее воз- 
можных значений числа классов в искомом нечетком с -разбиении Р”, где с, — наимень- 


шее, а с” — наибольшее из значений множества С. 
В [15] было предложено объединить ЕСМ-алгоритм с процедурой вычисления 
соответствующего показателя оптимальности числа нечетких классов Г.(Р) в иско- 


мом нечетком с -разбиении Р” в рамках одной процедуры, параметрами которой явля- 
ются значения наименее возможного с. и наиболее возможного с’ числа нечетких кла- 
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стеров в искомом нечетком с-разбиении Р”. Вместе с тем, предложенная в [15] про- 
цедура предъявляет достаточно высокие требования к оперативной памяти ПЭВМ, 
особенно при обработке совокупностей сравнительно большого объема, и в [7] была 
предложена модификация предложенной в [15] процедуры, получившая, от англоязычных 
терминов #177у с -теапз и сазег уаПацу, название ЕСМ-СУ-алгоритма. Рассмотренный 
в предпринятом исследовании метод построения множества наиболее возможного 


числа нечетких кластеров в искомом нечетком с-разбиении Р” позволяет модифи- 
цировать ЕСМ-СУ-алгоритм, и соответствующая модификация процедуры будет 
рассмотрена в процессе дальнейшего изложения. 


Построение множества значений возможного числа 
нечетких кластеров на основе нечетких чисел 


Для дальнейшего рассмотрения представляется необходимым напомнить опре- 
деления понятий нечеткой величины, нечеткого интервала и нечеткого числа. Если 
некоторое нечеткое множество Г определено на множестве действительных чисел, 
то есть представляет собой отображение 5% — [0/1], то Г именуется нечеткой величиной. 


Каждое значение х, е будет называться модальным значением нормальной нечет- 
кой величины У ‚если х, является элементом ядра Соге(Г) нечеткой величины Г ‚то 
есть и,(х,) =1. В случае, если саг4(Соге(у))= 1, то Г является унимодальной нечет- 
кой величиной, в случае же, когда саг(Соге(У))> 1, нечеткая величина Г называется 
мультимодальной [16]. 

Пусть Ё или К -— невозрастающая функция %* [0,1], такая, что [(0) = (0) =1 
и Ух, > 0, [(х,) < 1, Ух, <1, (х,)>0; [40] =0 либо имеет место Г(х,) > 0, Ух, и [(+0)=0. 
Тогда нечеткое множество Г называется нечетким интервалом ГА -типа с функцией 
принадлежности му (х,) , определяемой формулой 


<т, (7) 


где т называется нижним модальным значением 7, т — верхним модальным значе- 


нием, а параметры а > 0 и Ь>0 называются левым и правым коэффициентами нечет- 
кости соответственно. Таким образом, нечеткий интервал Г.К -типа может быть пред- 
ставлен в виде четверки параметров, что записывается в виде Г =(т‚т‚а,Ь) р. 


Если Г =(тт,а,Б),к является нечетким интервалом ГК-типа, то при условии 
совпадения нижнего и верхнего модальных значений т=т =т нечеткий интервал 


[К -типа У именуется нечетким числом [К -типа с функцией принадлежности, опре- 
деляемой в соответствии с выражением 


"я } при х,<т 


у(х,) = (8) 
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где т называется модальным значением нечеткого числа ГА-типа, символически 
записываемого в виде Г =(т‚а,Б) р, аа и Б — левым и правым индексами нечет- 


кости соответственно. 
При использовании оптимизационных методов нечеткой кластеризации число 


классов с в искомом нечетком с -разбиении Р” может определяться с помощью ре- 
зультатов разведочного анализа данных, к примеру, при анализе диаграммы 


рассеивания [2] или экспертных оценок. Пусть &, = {(&,, м в (&))|2=&, <п,\Уе= Ьй} = 


множество одноточечных нечетких множеств, выражающих мнения А экспертов о 
числе с кластеров в искомом нечетком с -разбиении Р”, так что предполагаемому е -м 
экспертом числу &, нечетких кластеров в Р” соответствует степень уверенности 
МЕ, (Е) = (0,|. Для каждого е=1.....К значение & может рассматриваться как мо- 
дальное значение нечеткого числа, и для каждого из строящихся нечетких чисел 
Г, =(т..а,,Б,) в, е=1...,К © модальными значениями т, = &, соответственно значе- 


ния у (1), е=1,...,К будут полагаться равными нулю в точках с=1 и с=п, так что 
и, (В=и, (п)=0, УееЦ,...,К}, левый коэффициент нечеткости будет определяться 


выражением а, = (Е, —1), а правый — Б, =(п- Е,), так что носителем 5ирр(У,) нечеткого 


числа И, е=|,....к будет открытый интервал (Е, -а,, Е +Ь,)=(1, п), а ш, (1) будет 


определяться выражением 


Е 
И - ‚Иры 152 


Ау, (= , > (9) 
К “| 1>Е 
г при 


где индекс е функций Ги К подчеркивает, что они могут выбираться для каждого 
нечеткого числа Г, =(т,,а,,6,) в, е=1....,К отдельно. Таким образом, если на коор- 


динатной плоскости по оси абсцисс откладывать число классов [е (1, п), а по оси ор- 
динат — значения функции принадлежности и, (Г) нечеткого числа Г,, то функция 
представления формы нечеткого числа У,, е=1,..., будет иметь вид кривой, дости- 
гающей максимума в точке с координатами /=&,, му. (1) =1. При построении нечетких 
чисел Г,, е=1....,К для того, чтобы каждое Г,, е=1,...,К описывалось бы непрерывной 
функцией принадлежности Ди. (1) ‚ подразумевается, что /е(Ьп)<\%. 

Далее следует построить нечеткую величину Г с непрерывной функцией при- 
надлежности д,([), для чего ко всем Г, =(т,,а,,Б,)», е=1.....К можно применить 


операцию объединения, причем в данном случае операция объединения понимается 
в широком смысле, то есть может быть осуществлена с помощью какой-либо 
выбранной исследователем $-нормы [17]. После построения У выбирается порог 
а = пт Ме. (Е), для которого строится И.) = Ми) (1))| Ми (Г) =, (1) > а} -— нечеткое 


множество уровня а с непрерывной функцией принадлежности Ми (Г). Носитель не- 
четкого множества И, в таком случае будет представлять собой интервал действи- 


тельных чисел $ирр(И „›)=[с", с"], и множество С = {с.,...,с”} возможных значений 


«Штучний 1нтелект» 32008 р 


Вятченин Д.А. 


числа классов в искомом нечетком с -разбиении Р” может быть определено как под- 
множество натуральных чисел, определяемое выражением 


{1 е, <", с}, (10) 
где [с] =с, — округление числа с' до ближайшего сверху целого из интервала [с",с"], 
[с"|=с* — округление числа с" до ближайшего снизу целого из интервала [с",с"], а 
7’ — множество положительных целых чисел. Для элементов с, множества допус- 


тимых значений числа нечетких кластеров С = {с,,...,с”} можно определить значения 


принадлежности ие(с,)= и У (1),1=с,,с, ЕС так что множество С будет представлять 


собой носитель нечеткого множества = {с,, Ик(с,)}, с, ЕС с дискретной функцией 
принадлежности. Значения функции принадлежности ие(с,), с, =С могут интерпре- 
тироваться как степени адекватности значений с, еС, так что обобщенный пока- 
затель оптимальности числа нечетких классов Г.(Р) может быть определен в виде 
Г. (Р)=И.(Р). Ив(сь) ‚ где с. ЕС — число кластеров в Р. 


Для построения множества одноточечных нечетких множеств С, е=1,....К можно 


воспользоваться предложенным в [18] эвристическим О-АЕС-ТС-алгоритмом нечет- 
кой кластеризации с выбором различных расстояний между нечеткими множествами. 
Результатом работы О-АЕС-ТС-алгоритма является распределение К*(Х) объектов 


исследуемой совокупности Х по априори неизвестному числу с,2<с<п нечетких 


@ -кластеров для некоторого вычисленного порога сходства а, позволяющего коли- 
чественно оценить наименьшую степень сходства объектов в нечетких а -кластерах 
распределения А*(Х). Полученное в результате работы О-АЕС-ТС-алгоритма при 


некотором выбранном расстоянии число с нечетких & -кластеров в А *(Х) может за- 
даваться как &,, а значение порога сходства х полученного распределения К *(Х) — 
как значение д (,). 


Схема модифицированного ЕСМ-СУ-алгоритма 


Построение нечеткого множества Е= {сьше(с,)}, с, ЕС с помощью О-АЕС-ТС- 
алгоритма либо его модификаций является этапом, предваряющим построение нечет- 
кого с-разбиения Р”, оптимального в смысле выбранного показателя Г, (Р). 

Общая схема предлагаемой модификации ЕСМ-СУ-алгоритма, основанной на 
вычислении У,(Р) для всех с, еС, и которую можно обозначить как (п)ЕСМ-СУ- 
алгоритм, выглядит следующим образом. 

1. Полагается с :=с. и с,:=с и значения числа классов с в искомом Р” 


упорядочиваются следующим образом: 2<с, <...<с,<...< мы © 


2. Полагается (:=1. 
3. Вычислять: 
3.1. с помощью ЕСМ -алгоритма вычисляется нечеткое с -разбиение Р(с,) 


исследуемой совокупности на с, классов; 
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3.2. вычисляется значение показателя У,(Р,)=Г,(Р,). Ик(с,) для полученного 


нечеткого с -разбиения Р(с,); 
3.3. производится проверка условия: 
если (<2, 
то осуществляется переход на шаг 5, 
иначе осуществляется переход на шаг 4. 
4. Производится проверка условия: 
если У, (Р‚) >И, (Р*), 
то осуществляется переход на шаг 5, 
иначе осуществляется переход на шаг 6. 
5. Полагается Р(с,)=Р*. 
6. Производится проверка условия: 
если (<р, 
то полагается (= (+1 и осуществляется переход на шаг 3, 
иначе нечеткое с -разбиение Р* является искомым результатом и алгоритм 
прекращает работу. 


Схема (т)ЕСМ-СУ-алгоритма построена, исходя из предположения, что в качест- 
ве Г.(Р) используется коэффициент разбиения (3), так что в случае, когда в качестве 


Г.(Р) используется другой показатель, для которого решение задачи определения 


оптимального числа классов отыскивается в виде ши, (Р)), с=2,....п-1, на шаге 4 
Е 


представленной выше схемы следует производить проверку условия Г, (Р,) <Г.(Р*); 
следует также указать, что (т)ЕСМ-СУ-алгоритм требует сохранения в оперативной 
памяти ПЭВМ только двух нечетких разбиений — вычисленного Р(с,) и текущего Р*, 
вместо множества П(с,,с”)= {Р(с,) [6 =1,..., р} возможных решений задачи классифи- 


кации. При и*(с,)=1Ус, ЕС, предложенная схема будет являть собой ЕСМ-СУ- 


алгоритм [7]. 


Иллюстративный пример 


Для проведения вычислительного эксперимента были выбраны изображенные 
на рис. 1 двумерные данные, предложенные в качестве тестовых К.Г. Луни [19], представ- 
ляющие собой совокупность 15 объектов Х = {51,...,х!5}. При проведении вычисли- 
тельных экспериментов исходные данные были пронормированы в соответствии с 
выражением 

ий (р 


тах & 
7 


и обработаны О-АЕС-ТС-алгоритмом с использованием относительного обобщен- 
ного расстояния Хемминга, относительного евклидова расстояния и относительной 
евклидовой нормы [18]. 


«Штучний 1нтелект» 32008 529 


Вятченин Д.А. 


НЫ 
0 1 2 Е 4 5 Б 7 Е эх 
Рисунок 1 = Данные ДЛЯ проведения вычислительного эксперимента 


При использовании относительного обобщенного расстояния Хемминга было 
получено распределение по & =2 при а, =0.8125, а при использовании относитель- 
ного евклидова расстояния и относительной евклидовой нормы было получено 
& =& =3 при а. =0.8065 и а; =0.9626 соответственно, что дает возможность пост- 
роить нечеткие числа И, = (т, а, В) в и И, =(т,,а>,Б.) к» с модальными значениями 
т =2, т, =3 и функциями принадлежности 


- 5 и =) 1</<2 
ыы 2 


ми (0 = 2+5 (12) 
а 5 5 ] 2<1<15 
и 
от ) 15/3 
(13) 
о ]} 35115 


1 


Рисунок 2 — Функция принадлежности нечеткой величины Г 
и значения принадлежностей нечеткого множества 
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Для построения Г в качестве $-нормы была выбрана операция объединения 


нечетких множеств, а значение х было выбрано как а=шта»,е=1,2,3 и составило 
е 
> =0.8065, так что множество С представляет собой совокупность С = {2,3,4,5,6} . 


На рис. 2 непрерывной кривой изображена функция принадлежности у (1) , а симво- 


лом о обозначены значения Ив(сь), с, ЕС нечеткого множества К. 


На рис. 3 изображено поведение для построенного обобщенного индекса 
разделимости Г,(Р) при обработке данных (т)ЕСМ-СУ-алгоритмом. 


Е (Р и. (Р) 
а ) | | 0.247 | 


Рисунок 3 — Значения И. (Р) (а) и Г,(Р) при Шве(с) =Ь Ус, ЕС (6) 


при обработке исходных данных (т)ЕСМ-СУ-алгоритмом 


На рис. 4 изображено поведение для построенного показателя ,.(Р). 


В обоих случаях оптимальным числом классов в искомом нечетком с -разбиении 
является с=3. Значения принадлежностей объектов совокупности Х={х,,...,х,5} 


первому нечеткому кластеру оптимального Р* изображены на рис. 4 символом о, 
второму — символом Д и третьему — символом м. 
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Рисунок 4 — Значения принадлежностей объектов классам оптимального Р* 


Анализ исходных данных, приведенных на рис. 1, позволяет визуально выде- 
лить три довольно четкие совокупности точек, и приведенные на рисунке 4 значения 
принадлежности объектов исследуемой совокупности трем классам нечеткого с -раз- 
биения, оптимального в смысле критерия (2), соответствуют относительной близости 
объектов прототипам нечетких кластеров. Кроме того, результаты проведенного 
эксперимента показывают корректность определенного в результате работы (п)ЕСМ- 
СУ-алгоритма числа нечетких кластеров в искомом Р*, оптимальном как в смысле 
индекса разделимости Г,(Р) , так и в смысле коэффициента разбиения Г, (Р). 


Заключение 


Результаты проведенного исследования наглядно демонстрируют, что аппарат 
нечетких чисел является эффективным средством для построения множества значе- 
ний наиболее возможного числа классов в искомом нечетком с -разбиении, а метод 
построения нечетких чисел для решения указанной задачи является достаточно 
простым и вполне объяснимым с содержательной точки зрения. Кроме того, предло- 
женный подход к построению нечетких чисел на основе задаваемых модальных 
значений является гибким в том смысле, что функция представления формы нечеткого 
числа, с одной стороны, а также $-норма для построения нечеткой величины -— с 
другой могут быть выбраны в зависимости от условий задачи. 

Так как ЕСМ-алгоритм, отыскивающий минимум функционала (2), представ- 
ляет собой параметрическое семейство по у при фиксированном числе кластеров с [2], 
и при увеличении значения у возрастает неопределенность классификации, что в свою 
очередь влияет на поведение показателей Г, (Р) ‚ то при больших значениях у иногда 
оказывается невозможным определить локальный экстремум показателя Г,(Р). Этим 
обстоятельством диктуется целесообразность использования обобщенного показателя 
Г.(Р) в предложенном (п)ЕСМ-СУ-алгоритме, анализ данных с помощью которого 
позволяет в значительной мере снизить требования к оперативной памяти ПЭВМ при 
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обработке больших массивов данных. В свою очередь, главным достоинством предло- 
женного двухэтапного подхода к решению задачи нечеткой кластеризации, заключаю- 
щегося в совместном использовании О-АЕС-ТС-алгоритма и (т)ЕСМ-СУ-алгоритма, 
является возможность обработки данных в полностью автоматическом режиме. 
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Д.А. Вятчешн 

Використання неч1тких чисел задля обгрунтування кластер в у методах нечтко! кластеризаци 
У статт! розглядаеться метод визначення оптимального числа кластер1в у нечткому с-розбитть, 
заснований на побудов! 1нтервалу значень найблльш можливого числа класв з використанням 
нечиких чисел. Запропонована модифжашя ЕСМ-СУ-алгоритму 1 наводиться результат обчислювального 
експерименту. 
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